Time domain radiation and absorption by subwavelength sources 
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Radiation by elementary sources is a basic problem in wave physics. We show that the time- 
domain energy flux radiated from electromagnetic and acoustic subwalength sources exhibits re- 
markable features. In particular, a subtle trade-off between source emission and absorption underlies 
the mechanism of radiation. This behavior should be observed for any kind of classical waves, thus 
having broad potential implications. We discuss the implication for subwavelength focusing by time 
reversal with active sources. 
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Any textbook on wave physics or field theory contains 
a chapter on radiation by elementary sources in homo- 
geneous media ,l,j. Frequency-domain analyses of the 
radiated fields and the associated energy fluxes are the 
most widespread. In these approaches, the far-field en- 
ergy flux is usually defined as the contribution that sur- 
vives time averaging, corresponding to power that contin- 
uously leaks away from the source. Conversely, radiated 
near fields generate oscillating terms in the energy flux, 
that are discarded in the time-averaging process. Time- 
domain expressions of radiated fields are also common in 
the context of electromagnetic radiation J3j5| , including 
the optical regime [y], and in acoustics [3-Q- Neverthe- 
less, time-domain expression of the energy flux have been 
given much less consideration. In the case of electromag- 
netic waves, the time-domain energy flux radiated from 
an electric dipole at rest may be found in some textbooks 
(see Ref. 3 for instance). Its expression is also at the 
core of interesting studies of the time decay of classical 
oscillating dipoles [Ifl, \lM ; but that do not describe the 
full contribution of the near-field terms that is discussed 
in the present study. In most textbooks, the discussion is 
limited to harmonic oscillations and time averages, since 
the focus is usually on far-field radiation 0, Q . In acous- 
tics, although time-domain expressions of the fields radi- 
ated by monopole or dipole sources are widespread [7H9| , 
we are not aware of any discussion of the time-domain en- 
ergy flux, and in particular of its near-field and far-field 
components. 

In this Letter, we revisit the basic problem of radiation 
by elementary subwavelength sources, from the point of 
view of emission and absorption of energy in the time do- 
main. Considering time-domain expressions of the energy 
flux for the acoustic monopole and the electromagnetic 
dipole, and analyzing carefully the energy balance, we 
show that there is a subtle trade-off between emission of 
energy and subsequent reabsorption by the source, the 
difference between emission and reabsorption giving the 
amount of energy that is irreversibly radiated to the far 
field. This result reveals some important features of the 
dynamic interchange of energy between a subwavelength 



source and a wavefield, that have not been discussed so 
far, to the best of our knowledge. It also suggests a 
novel point of view on near-field radiation. Since the 
conclusions hold for both acoustic and electromagnetic 
waves (with striking similarities), they underline a be- 
havior that should be found with any kind of classical 
waves, thus having broad implications. We illustrate an 
implication in the context of subwavelen gth focusing us- 
ing time reversal with active sources jlj, [l3| • 

The propagation of electromagnetic waves generated 
by a spatially localized source in an otherwise homo- 
geneous medium is described by the following equa- 
tion 0, i] 
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(r,i) + Vx VxE(r,i) = Se™(r,i) (1) 



where E(r,t) is the electric field at point r and time t, 
and c is the speed of light in the medium. The source 
term Sem(r, i) is often written in the form Semi^,t) = 
—fio (9/9i)j(r, t), where j(r, t) is the electric current den- 
sity and /io the vacuum magnetic permeability. The 
electromagnetic energy current is given by the Poynt- 
ing vector n(r, i) — E(r,t) x H(r,t), where E(r, i) is 
the retarded solution of Eq. ([IJ and H(r,t) the associ- 
ated magnetic field. The energy flux (j)em{R,t) across a 
sphere with radius R centered at the origin is 4>em {R, t) ~ 
/sphere n(r,i) ' udV, where u = r/|r|. 

For acoustic waves in the linear regime, the acoustic 
pressure field p(r, i) generated by a spatially localized 
source in a homogeneous medium obeys [3, Q : 
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where Cg is the acoustic velocity in the medium and 
5(jc(r,i) the source term. The acoustic energy current is 
q(r, t) — p(r, t)v(r, i), p(r, t) being the retarded acoustic 
pressure field solution of Eq. ([2]) and v(r,i) the associ- 
ated acoustic velocity field. The energy flux follows from 

'/'ac(i?,0 = /.phcrcq(l-,i)-udV 

In this Letter we study the radiation produced by 
sources of size much smaller than the characteristic 
length scale of the wavefield, that will be denoted by 
"subwavelength sources" . In the case of electromagnetic 



waves, we use a point electric dipole model, with dipole 
moment p(i) = f{t)po, f{t) being the dimensionless 
time-domain amplitude and po a time-independent vec- 
tor accounting for the source polarization. This model 
describes, e.g., a dipole moment p(i) = (7eL(t) corre- 
sponding to an oscillating charge qe with oscillation am- 
plitude L(i) much smaller than all other relevant charac- 
teristic lengths 3]. For a dipole centered at r = 0, the 
electromagnetic source term reads: 



bem(r, t) 






(3) 



where ^(r) is the three-dimensional Dirac delta function. 
In the case of acoustic waves, we use a point mass source 
model describing a radially oscillating sphere with radius 
a{t) = oo -\- ^{t), in the limit of of vanishingly small ra- 
dius Q. For a source centered at r = 0, the acoustic 
source term reads: 



Sac(j, t) = Po So 
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where po is the mass density of the unperturbed homo- 
geneous medium and sq = AnaQ. For the sake of formal 
similarity with the electromagnetic case, we will write 
^{t) = f{t)(,o with ^0 a- time-independent length driving 
the acoustic source strength. 

The time-domain solutions of Eqs. ([IJ and ([2]) with the 
source terms given by Eqs. ^ and (JH) can be found in 
textbooks on electromagnetic and acoustic waves propa- 
gation @t3;H13- From the field expressions, the energy 
flux across a sphere with radius R can be deduced af- 
ter tedious but straightforward algebra. In the case of 
electromagnetic waves, one obtains: 
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For acoustic waves, the explicit calculation of the energy 
flux leads to: 



AR.t) = 
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di 
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.(6) 



In Eqs. (O and ^ all terms within square brackets [...] 
denote retarded values, and have to be evaluated at time 
t — R/c (electromagnetic waves) or i — R/cg (acoustic 
waves). Although their derivation is a rather simple ex- 
ercise, we will see that these expressions bring to light 
fundamental aspects of the mechanism of radiation by 
subwavelength sources that have not been discussed so 
far. 

^From a qualitative point of view, the structure of 
Eqs. ([U and ([6|) deserves several comments. The far- 
field limit, obtained for R —^ oo, leads in both cases to 
an energy flux proportional to the square of the second 



derivative of the source amplitude, in agreement with a 
well-established result in classical wave theory [1| . For a 
monochromatic source oscillating at a frequency w, with 
f{t) — sm{Lut), this far-field term is the only one that 
survives a time-averaging of Eqs. (O and ([6]). The far- 
field behavior is extensively discussed in textbooks, both 
for monochromatic and pulse sources. Nevertheless the 
time-domain electromagnetic and acoustic energy fiuxes 
contain additional near-field terms whose amplitude de- 
pend on the distance R to the source. The first near-field 
term scales as R~^ and is identical in Eqs. dSJ and ([6|), 
except for a factor of two, while additional terms scal- 
ing as i?~2 and R~^ appear only in the expression for 
the electromagnetic case. These near-field contributions 
exhibit remarkable properties that induce specific behav- 
iors of the time-domain energy fiux. A first result is that 
the time-dependent amplitudes of the near-field terms in 
Eqs. ([U and ^ read as first-order derivatives of func- 
tions that are positive (squares) and that recover their 
initial values after a finite time interval (the pulse dura- 
tion, or the period for monochromatic excitation). As a 
result, these amplitudes necessarily change sign during 
their time evolution, meaning that the near-field terms 
lead alternatively to outgoing or incoming contributions 
to the energy fiux. Conversely, the far-field term only 
contributes to an outgoing energy fiux. While this seems 
to be a commonly accepted result in the harmonic regime 
(for electromagnetic waves, it is known that the Foynting 
vector in the near field changes sign during one cycle of 
oscillation), the above result precisely demonstrates that 
the change of sign in the near-field energy flux also ex- 
ists for a pulsed source with finite duration (i.e., with an 
amplitude starting from zero and vanishing after a finite 
time). 

In order to study the behavior of the time-domain en- 
ergy flux on a quantitative basis, we need to specify the 
source amplitude function f{t). In the present work, we 
consider pulses with two requirements. First, f{t) has to 
be of strictly finite duration (denoted as T in the follow- 
ing), in order to define exactly a pulse onset {t = 0) 
and a pulse end (i = T). Second, f{t) and its time 
derivatives have to vanish continuously to zero at i = 
and i = T, in order to avoid temporal singularities in 
the energy fiux, as made clear from Eqs. ([S]) and ^. 
A broadband pulse matching these two requirement is 
for instance f{t) = exp[2T^/{t{t - T))] for t &]0,T[ and 
f{t) = otherwise. The temporal shape of the source 
amplitude f{t) and the shape of the associated far field 
amplitude are shown in Fig. [T] For such a broadband 
pulse, the period is on the order of the duration. More 
precisely, for the function f(t) given above, the period is 
close to half the duration, and the corresponding wave- 
length is A = cT/2 (see Fig.[T]). In this expression and in 
the following, for sake of brevity and since the velocities 
play the same role in the electromagnetic and acoustic 
cases, both c and Cs are referred to as c. 

The knowledge of f{t) and its derivatives allows us to 
plot the time evolution of (j)em{R, t) and 4>ac{R, t) for dif- 
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FIG. 1: Time evolution of the source amplitude f{t) (left) 
and of its second derivative d^/(t)/di^ (right). The latter 
represents the time-dependence of the far-field amplitude for 
both electromagnetic and acoustic waves. 



ferent observation distances, covering the near-field, the 
intermediate and the far-field regimes. We show in Fig. [2] 
the time evolution of the energy flux in the electromag- 
netic (top) and acoustic (bottom) situations, and for four 
different distances. In the far field (i? ^ A), the energy 
flux is always positive and describes the radiated energy 
flowing irreversibly from the source. In the near field 
(i? ^ A), a completely different behavior is observed. 
The energy flux oscillates, and takes negative values on 
some time intervals. This means that part of the energy 
that has flowed outside the sphere of radius i? at a given 
time flows back into the sphere at subsequent times. 

At this stage, conservation of energy states that a neg- 
ative energy flux corresponds to an increase of energy 
stored inside the sphere with radius R, or to reabsorp- 
tion into the source (or both) . In order to quantitatively 
settle this point, we introduce Ux{R, t) defined as the en- 
ergy stored outside the sphere with radius R at time t in 
the electromagnetic or acoustic field (the subscript "x" 
stands for em or ac). It reads: 



Ux{R,t)^ / 0,(i?,t')dt' 



(7) 



The time evolution of Uem {R, t) is shown in Fig. [3] for 
the same distance regimes as in Fig. [5] Although not 
shown for the sake of brevity, the same behavior is ob- 
served for acoustic waves. As expected from the changes 
in sign of the energy flux, we see that Uem {R, t) is not a 
monotonic function of time except in the far field. This 
non-monotonic behavior of the time evolution of the en- 
ergy stored in the field can be characterized by splitting 
Ux{R, t) into U.x{R, t) = U^ + AC/^(i?, t). The first term 
C/^ = /p°° 03;(i?, i) dt corresponds to the overall time- 
averaged energy eventually radiated irreversibly through 
the sphere of radius R to the far field, and is independent 
of R. The second term describes the time variations of 
the energy stored in the field beyond the distance R, and 
either increases or decreases Ux {R, t) with respect to the 
asymptotic value U^ . This dynamic behavior is fully 
described by the curves in Fig. [31 One clearly sees that 
at some time range, for i? <C A, the energy stored outside 



the sphere with radius R exceeds the final energy that 
remains in the field after the source has been turned off 
[t > T). This proves that part of the energy of the field 
has been reabsorbed by the source, which constitute the 
main result of this work. This result, derived here using 
a pulse of finite duration and finite energy, remains valid 
for monochromatic and quasi-monochromatic waves. It 
shows without ambiguity that a negative energy flux ob- 
served in the near field corresponds to reabsorption by 
the source. 

This conclusion puts forward new features of the near 
field. Although it is known that on average, near-field 
terms correspond to non-radiative energy [7|, Il4l . Il5| , our 
work shows that this non-radiative energy is dynami- 
cally exchanged between the field and the source, at 
the time scale of the main oscillation. This subtle dy- 
namic process is hidden in the first-place when compu- 
tations are restricted to time-averaged values. We also 
stress that a time-domain analysis reveals behaviors that 
cannot be seen in the frequency domain. For example, 
in near-field optics or acoustics, it is often stated that 
some information is lost in the far field due to the loss 
of non-radiative components that remain spatially local- 
ized close to the sources (in the near field zone). With a 
non-stationary source, one could question what happens 
after the source has been turned off. Is the field finally 
radiated into the far field, and if so, where is the loss of 
information? Our work provides an unexpected answer: 
in the near field, some energy is constantly dynamically 
exchanged between the field and the source, and even- 
tually most of it is absorbed by the source while only 
a small part is radiated into the far field. The discus- 
sion has been limited in this study to a subwavelength 
source emitting in a homogeneous medium, so that only 
near fields produced by the source itself have been con- 
sidered. A more general analysis including near fields 
produced by scattering from subwavelength objects (sec- 
ondary sources) should also reveal interesting dynamic 
behaviors. In particular, understanding, in the time do- 
main, the concept of non-radiative components that ap- 
pear in the frequency-domain angular spectrum decom- 
position of scattered fields [ia| would be another step 
forward. This is left for future work. 

It is also interesting to have a look at the distance 
dependence in the near field of the maximum value of the 
energy stored in the field AC/™''^(i?) = jaayi{MJx{R, t)}. 
Conserving only the dominant terms as i? — > in Eqs. ([5]) 
and dH), it is easy to show that AU^''{R) ~ R-^ and 
AC/™°^(i?) ~ R~^. Therefore for a quasi point source 
model, the energy transiently stored in the field becomes 
arbitrarily large at short distance. In practice, the energy 
must be limited somehow by the limitations on the source 
model itself. Another peculiar behavior, observable only 
with broadband pulses of strictly finite duration, is that 
the energy fiux exhibits a slight sign inversion even at 
times t > T, i.e., after the source has become inactive 
(see the insets in Fig. [5] and [3]). This sign inversion does 
therefore not correspond to reabsorption in the source 
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FIG. 2: Time evolution of the electromagnetic and acoustic energy flux (j)em.{R, t) (top row) and cf>ac{R, t) (bottom row) for four 
different distance regimes. Far-field regime R^ \, limit of the source free regime R — 2X (2A = cT in this particular case), 
near-field regimes R < X and R <^ X. For R = 2A, the insets show the sign inversion of the energy fiux. 
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FIG. 3: Time evolution of the electromagnetic energy (7em(-R, f) stored outside the sphere of radius R at time t, for the same 
distance regimes as in Fig. (2] The inset shows the dip due to the sign inversion of the energy flux. 



in this case, but to a small part of the energy flowing 
back and forth through the sphere of radius R. This 
"anomaly" becomes insignificant (although non strictly 
zero) in the far field since it is due to the contribution of 
terms in the energy flux that decay as R^^ or faster. 

To our knowledge, the near-field contributions in the 
time-domain energy flux have been first discussed in op- 
tics by Mandel [10], in the context of the decay rate of 
a classical electric dipole in vacuum. The discussion was 
constrained by the fact that for a freely decaying atomic 
dipole, "the total field energy could not exceed the max- 
imum amount of energy of the dipole that ultimately 
emerges by radiation" flO|. A major difference with the 
present work is that Mandel's approach considered the 
time variation of the envelope of the emitted wavefield, 
but terms varying in time at the scale of the optical pe- 
riod were discarded. In a more recent study, Schantz 
considered the time variations of the energy flux emitted 
by a decaying electric dipole keeping all time-dependent 
terms, with an initial condition corresponding to an elec- 
trostatic dipole [11|. He concluded that the eventually 



radiated energy had to correspond to electrostatic energy 
initially stored in the far field. Although it is out of the 
scope of this Letter to further discuss this unexpected 
and interesting result, we point out that the situation 
studied by Schantz is very different from that considered 
here. Indeed, we considered as a fundamental assump- 
tion the case of a medium initially free of energy, with a 
source amplitude starting exactly from zero, and vanish- 
ing rigorously after a finite time (as opposed to an initial 
non-zero static field). 

The results presented in this Letter were derived in 
the case of electromagnetic and acoustic radiation, but 
they certainly underline general behaviors that should 
be found for any classical subwavelength sources. There- 
fore, the peculiar dynamics of the energy exchange be- 
tween a subwavelength source and the radiated field has 
potentially broad implications. Here, we discuss one im- 
portant consequence in the context of subwavelength fo- 
cusing by time reversal. Experimental realizations of 
time-reversed wavefields have been demonstrated both in 
acoustics and electromagnetism, by use of close 2D or 3D 



cavities [Ij, [13, [3 ■ When the field emitted by a point- 
hke source is time-reversed in the source-free medium, 
refocusing is hmited by diffraction [l7|. However, when 
both the wavcficld and the source are time reversed, per- 
fect refocusing can be obtained [ij, [l6| . Accordingly, ex- 
periments in acoustics have demonstrated subwavelength 
refocusing with an active time- reversed source |13| . the 
focal spot size being limited only by the finite size of the 
source itself. Intuitively, the role of the time-reversed 
source is seen as that of a sink, i.e., of an absorber of 
the incoming time-reversed wave. Our work shows that 
the role of the time-reversed source is more subtle, and 
that it necessarily involves both absorption and emission 
of energy. Indeed, the time-domain evolution of the field 
energy in a perfect time-reversal experiment (with re- 
versed field and source) is directly given by the curves in 
Fig. 3 read backwards. Therefore, the energy in the field 
is transiently larger than the energy carried by the time- 
reversed wavefield, so that in some time range, the sink 
actually behaves as a source. The time-reversed source 
is both an absorber and an emitter. The term "sink" 
therefore only makes sense when one considers the over- 
all energy balance, obtained after time integration. Our 
work has two important consequences for practical ex- 
periments. First, the focusing performances cannot be 
discussed without considering the energy point of view, 
in particular because for a sink of vanishingly small size, 
the transient energy that has to be stored in the field be- 



comes arbitrarily large. Second, perfect subwavelength 
refocusing (i.e., without energy scattered away from the 
focal spot) cannot be achieved by use of a passive sub- 
wavelength absorber, as efficient as it may be, since the 
dynamic exchange of energy is a necessary condition for 
a localized absorption of the full energy of the wavefield. 

In summary, from the study of time-domain expres- 
sion of the energy flux radiated by pulsed electromag- 
netic and acoustic elementary sources, we have shown 
that the non-radiative energy predominant in the near- 
field is dynamically exchanged between the source and 
the field. We have discussed implications for subwave- 
length focusing and imaging. Since the results hold for 
both electromagnetic and acoustic waves, we believe that 
they underly a universal process of radiation by any kind 
of subwavelength sources, although demonstrated here 
only for the acoustic monopole and the electromagnetic 
dipole. In the case of electromagnetic waves emitted by 
a single classical dipole emitter, a giant transient storage 
of electromagnetic energy is necessary in order to radiate 
a (much smaller part) in the far field. It would be in- 
teresting to clarify the way quantum theory handles this 
point in the computation of spontaneous emission by a 
single atom. 
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